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Abstract. A very simple and short proof of the polynomial matrix spectral factorization 
theorem (on the unit circle as well as on the real line) is presented, which relies on elementary 
complex analysis and linear algebra. 
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I. Introduction 

In this paper, we present an elementary proof of the polynomial matrix spectral 
factorization theorem: 

Theorem 1. Let 

N 

S(Z) = 

n=-N 

be an m x m matrix function (C n G <C mxm are matrix coefficients) which is positive 
definite almost everywhere on T := {z G C : \z\ = 1}. Then it admits a factorization 

S(z) = S + (z)S-(z), z G C\{0}, 

where S + (z) = Xl^lo^™ 2 ™ is an m y. m polynomial matrix function which is nonsin- 
gular inside T, det S + (z) ^ when \z\ < 1, and S~(z) = S + (l/z) T = YlH=o^*n z ~ n zs 
its adjoint, A* = A n T , n = 0,1, . . . , N , (respectively, S~ is analytic and nonsingular 
outside T). S + is unique up to a constant right unitary multiplier. 

In the scalar case, m — 1, the above result is known as the Fejer-Riesz lemma and 
can be easily proved by considering the zeroes of S(z). 

The matrix spectral factorization S(t) = S + (t) (i)) , \t\ = 1, was first established 
by Wiener [8] in a general case for any integrable matrix-valued function S(t), \t\ = 1, 
with an integrable logarithm of the determinant, logdetS'(t) G £i(T). In this case, 
the spectral factor S + belongs to the Hardy space H 2 . Wiener proved this theorem 
by using a linear prediction theory of multi-dimensional stochastic processes. A little 
bit later, by using the same methods, Rosenblatt [7] showed that S + is a polynomial 
whenever S is a Laurent polynomial. Since then, many different simplified proofs 
of the matrix spectral factorization theorem have appeared in the literature, see for 
example [9], [2], [1], [4]. A very short proof of this result is given in [3], however it 
uses Wiener's general spectral factorization existence theorem and some facts from the 
theory of the Hardy spaces. 

The presented proof relies only on elementary complex analysis and linear algebra. 
This proof is constructive, which makes it possible to compute the spectral factor 
approximately at least in the case of low dimensional matrices (see [5] for a new 
reliable computational algorithm of general matrix spectral factorization). The same 
pattern can be also used for proving in a straightforward manner the polynomial matrix 
spectral factorization theorem on the real line: 
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Theorem 2. Let 

2N 

S(z) = J2C n z n 

n=0 

be an m x m matrix function (C n G C mxm are matrix coefficients) which is positive 
definite almost everywhere on the real line R. Then it admits a factorization 

S(z) = S + {z)S-{z), zeC, 

where S + (z) = J2n=oAnZ n is an m x m polynomial matrix function which is non- 
singular in the open upper half plane, det S + (z) ^ when Ivciz > 0, and S~(z) = 
S + (z) T = J2n=o J ^n 2;n i s tt s adjoint, A* n = A n T , n = 0, 1, . . . , N, (respectively, S~ is 
nonsingular in the open lower half plane) . S + is unique up to a constant right unitary 
multiplier. 

Various practical applications of spectral factorization in linear systems are widely 
recognized (see, e.g. [6]), where the problem naturally arises in either of the two 
different forms commonly called discrete and continuous. Mathematically both forms 
are equivalent under a conformal mapping of the upper half plane into the unit disk. 

2. Notation 

For a G C, a* = a denotes its conjugate, and for a matrix A, A* = A T denotes its 
Hermitian conjugate. 

Let T + = {z G C : \z\ < 1}, T_ = {z G C : \z\ > l}U{oo}, R + = {z G C : Imz > 0}, 
and R_ = {z G C : Imz < 0}. 

Lf(T) (Lj~(T)) stands for the class of integrable functions on T whose Fourier 
coefficients with negative (positive) indices equal to zero. 

Let TZ mxrn be the ring of rational mxm matrix functions defined in the complex 
plane. For / G 7^ mxm - ; the adjoint matrix function / is defined by f(z) = f (l/z) T 
in the discrete case and by f(z) — f (~z) T in the continuous case. Since / is uniquely 
determined by its values on T (on R), and f(z) = (f(z))* for z G T (for z G R), usual 
relations for adjoint matrix functions, like fg(z) = g(z)f(z) and f~ l (z) = f~ l (z), etc., 
are valid. Obviously, if f(e ie ) G L+(T), then f(e id ) G L{ (T). Whenever S + G TZ mxm 
is determined, S~ always denotes its adjoint. 

U G TZ mxm is called paraunitary if 

U(z)U(z) = I m 

where I m stands for the m-dimensional unit matrix. Note that U (z) is a usual unitary 
matrix on the boundary, i.e. 

(1) U(z)U(z)* =I m , zeT (z G R). 

We say that a matrix function is analytic in a domain if the entries of the matrix 
are analytic in the domain. 
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3. An elementary proofs of Theorems 1 and 2 

Proof of Theorem 1. By Gauss elimination on the matrix S(z) and the Fejer-Riesz 
lemma, a factorization 

(2) S(z) = S (z)S (z) 

can be easily achieved with Sq(z) G TZ mxm , (Namely, if A = (o^Oi^Tm * s a P os i" 
tive definite matrix and B = (&ij)ij=Tm * s the unique positive definite matrix such 
that A = BB*, then the entries of B can be recursively determined by the for- 
mulas fen = y/cfa, b kl = Ofci/^/oIT*, fc = 2,3, ...,m, 6 nn = - Y^jZlKjKji 

bkn = {akn - Y^jZi hj^j /Kn> n = 2 , 3, . . . , to, = n + 1, . . . , to. These formulas re- 
main valid for rational matrix functions as well. We need only to assume that y/a and 
a* are the scalar spectral factor of a and the adjoint of a, respectively.) 

If Sij is the ijth entry of So with a pole at a G T + , then we can multiply S 
by the paraunitary matrix function U(z) = diag[l, . . . , u(z), . . . , 1], where u(z) = 
(z—a)/(l—az) is the jjth entry oiU(z), so that the ijth entry ofthe product S (z)U(z) 
will not have a pole at a any longer keeping the factorization (2): (SqU)(z) SqU(z) = 
So(z)Sq(z) = S(z). In the same way, one can remove every pole of the entries of So at 
points in T + . Thus S can be represented as a product 

(3) S(z) = S+(z)S (z), 

where Sq G TZ mxm is analytic in T + 

Now, it might happen so that Sq is not nonsingular everywhere on T + . If \a\ < 1 
and detS' f (a) = 0, then there exists an to x to unitary matrix U such that the 
product Sq(o)U has all 0's in the first column. Hence a is a zero of every en- 
try of the first column of the matrix function Sq(z)U and the product Si(z) = 
Sq(z)U di&g[u(z), 1, . . . , 1], where u(z) = (1 — az)/(z — a), remains analytic inside 
T. While the factorization (3) remains true replacing Sq and Sq by S± and Sf , re- 
spectively, the determinant of S^ will have less zeros in T + than the determinant of 
Sq. Thus, continuing this process if necessary, we can remove any singularities in T + 
and get the factorization 

(4) S(z)=S + (z)S-(z), 

where S + G TZ mxm is analytic and nonsingular in T + . 

Now let us show that S + is in fact a polynomial matrix function of order N. 
S + is free of poles on T since S + (z)(S + (z))* = S(z) for z G T, and Z - N S + {z) = 

z~ N S(z)(S~(z)) 1 is analytic in T_. Consequently S + is analytic in C and z~ N S + (z) 
is analytic in T_ which implies that S + is a polynomial of order N. 

The proof of the uniqueness of S + is standard and it is given only for the sake of 
completeness. Namely if S(z) = S± (z)Si (z) = S'^"(^)S'^"(^) are two spectral factoriza- 
tions of S(z), then (S^iz)) 1 S±(z) is an analytic and nonsingular in T + paraunitary 
matrix function, which is free of poles and singularities on T because of (1). There- 
fore (^(e if? )) _1 1 S 1 + (e ie ) G Lf(T) n L{ (T), which implies that it is a constant matrix 
function. 
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Proof of Theorem 2. This proof can be carried out directly by the same steps as 
in the discrete case. The factorization (2) can be performed exactly in the same way 
as in the proof of Theorem 1 assuming under y/a a scalar spectral factor of rational 
function a which is positive almost everywhere on R. Elimination of poles of Sq (see 
(2)) and singularities of Sq (see (3)) in R + can be made by using the paraunitary 
factors u(z) = and u(z) = |5f > respectively. When we have the factorization (4) 
where S + G TZ mxm is analytic and nonsingular in R + , we can prove that S + is in 
fact a polynomial as follows: S + is free of poles on R since S + (z) (S + (z))* = S(z) for 
z G R (as in the discrete case), and S + = S(S')^ 1 is analytic in R_. Consequently 
S + G TZ mxm is analytic in C and hence polynomial. The order of S + is N since if A 
ia a nonzero matrix coefficient of the highest order of z in S + , then A A* ^ is the 
matrix coefficient of the highest order of z in S + S~ . 

The problem of uniqueness of S + can be reduced to the discrete case by using the 
linear fractional transformation z —> jz§ which maps T + to R + . 
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